Biophysical modeling of hemodynamic-based neuroimaging techniques

Thesis (Ph. D. in Medical Engineering and Medical Physics)--Harvard-MIT Program in Health Sciences and Technology, 2013.Cataloged from PDF version of thesis.Includes bibliographical references (pages 163-182).Two different hemodynamic-based neuroimaging techniques were studied in this work. Near-Infrared Spectroscopy (NIRS) is a promising technique to measure cerebral hemodynamics in a clinical setting due to its potential for continuous monitoring. However, the presence of strong systemic interference in the signal significantly limits our ability to recover the hemodynamic response without averaging tens of trials. Developing a new methodology to clean the NIRS signal from systemic interference and isolate the cortical signal would therefore significantly increase our ability to recover the hemodynamic response opening the door for clinical NIRS studies such as epilepsy. Toward this goal, a new method based on multi-distance measurements and state-space modeling was developed and further optimized to remove systemic physiological oscillations contaminating the NIRS signal. Furthermore, the cortical and pial contributions to the NIRS signal were quantified using a new multimodal regression analysis. Functional Magnetic Resonance Imaging (fMRI) based on the Blood Oxygenation Level Dependent (BOLD) response has become the method of choice for exploring brain function, and yet the physiological basis of this technique is still poorly understood. Despite the effort, a detailed and validated model relating the signal measured to the physiological changes occurring in the cortical tissue is still lacking. Modeling the BOLD signal is challenging because of the difficulty to take into account the complex morphology of the cortical microvasculature, the distribution of oxygen in those microvessels and its dynamics during neuronal activation. Here, we overcome this difficulty by performing Monte Carlo simulations over real microvascular networks and oxygen distributions measured in vivo on rodents, at rest and during forepaw stimulation, using two-photon microscopy. Our model reveals for the first time the specific contribution of individual vascular compartment to the BOLD signal, for different field strengths and different cortical orientations. Our model makes a new prediction: the amplitude of the BOLD signal produced by a given physiological change during neuronal activation depends on the spatial orientation of the cortical region in the MRI scanner. This occurs because veins are preferentially oriented either perpendicular or parallel to the cortical surface in the gray matter.by Louis Gagnon.Ph.D.in Medical Engineering and Medical Physic

Dynamic state-space estimation of the hemodynamic response with Near-Infrared Spectroscopy

Thesis (S.M.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, 2011.Cataloged from PDF version of thesis.Includes bibliographical references (p. 68-74).Near-Infrared Spectroscopy (NIRS) allows the recovery of the hemodynamic response associated with evoked brain activity. The signal is contaminated with systemic physiological interference which occurs in the superficial layers of the head as well as in the brain tissue. The back-reflection geometry of the measurement makes the DOI signal strongly contaminated by systemic interference occurring in the superficial layers. A recent development has been the use of signals from small source-detector separation (1 cm) optodes as regressors. Since those additional measurements are mainly sensitive to superficial layers in adult humans, they help in removing the systemic interference present in longer separation measurements (3 cm). Encouraged by those findings, we developed a dynamic estimation procedure to remove global interference using small optode separations and to estimate simultaneously the hemodynamic response. The algorithm was tested by recovering a simulated synthetic hemodynamic response added over baseline DOI data acquired from 6 human subjects at rest. The performance of the algorithm was quantified by the Pearson R2 coefficient and the mean square error (MSE) between the recovered and the simulated hemodynamic responses. Our dynamic estimator was also compared with a static estimator and the traditional adaptive filtering method. We observed a significant improvement (two-tailed paired t-test, p < 0.05) in both HbO and HbR recovery using our Kalman filter dynamic estimator compared to the traditional adaptive filter, the static estimator and the standard GLM technique. We then show that the systemic interference occurring in the superficial layers of the human head is inhomogeneous across the surface of the scalp. As a result, the improvement obtained by using a short separation optode decreases as the relative distance between the short and the long measurement is increased. NIRS data was acquired on 6 human subjects both at rest and during a motor task consisting of finger tapping. The effect of distance between the short and the long channel was first quantified by recovering a synthetic hemodynamic response added over the resting-state data. The effect was also observed in the functional data collected during the finger tapping task. Together, these results suggest that the short separation measurement must be located as close as 1.5 cm from the standard NIRS channel in order to provide an improvement which is of practical use. In this case, the improvement in Contrast-to-Noise Ratio (CNR) compared to a standard GLM procedure without using any small separation optode reached 50% for HbO and 100% for HbR. Using small separations located farther than 2 cm away resulted in mild or negligible improvements only.by Louis Gagnon.S.M